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Abstract 

We present a new model independent amplitude analysis of reaction ir~p —* 
7r~7r + n measured at CERN at 17.2 GeV/c on polarized target using a Monte 
Carlo method in an extended range of dipion mass 580-1080 MeV. The two 
solutions for the measured moduli \S\ 2 T, and \ S\ 2 T, of the two S-wave transver- 
sity amplitudes S and S show resonant behaviour below 800 MeV corre- 
sponding to a scalar resonance a (770). Simultaneous fits to |5| 2 £ and \S\ 2 T, 
including /q(980) resonance give solution average m a = 778 ± 16 MeV and 
r CT = 142 ±33 MeV. The CERN data on polarized target supplemented by the 
assumption of analyticity of production amplitudes in dipion mass allow to 
determine the helicity amplitudes So (a% exchange) and S± (tt exchange) from 
the fitted transversity amplitudes. The sign ambiguity in <r(770) contribution 
leads to two solutions. In the "down" solution cr(770) is suppressed in the 
helicity flip amplitude Si and thus also in ir~Tr + — > tt^tt + scattering. Most 
contribution of <r(770) to pion production is in the nonflip amplitude So, or 
in TT~af — > 7r~7r + scattering. In the "up" solution the situation is reversed. 
The "up" solution is excluded by unitarity in tttt scattering. The "down" 
solution, and thus the evidence for <r(770), is in agreement with unitarity in 
both 7r~7r + — > 7r~7r + and ir~al — > -k~tt + scattering. There are four "down" 
solutions (1,1), (2,1), (1,2) and (2,2). The cr(770) resonance manifests itself 
as a broad resonant structure at ~ 720 MeV in the flip amplitude |Si| 2 in 
solutions (1, 1) and (2, 1). The contribution of cr(770) to |Si| 2 is small in solu- 
tions (1,2) and (2,2). The metamorphosis of a narrow <r(770) in the nonflip 
amplitude So into a broad resonant structure in the flip amplitude Si is a new 
phenomenon related to breaking of scale and chiral symmetry in QCD. 
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I. INTRODUCTION 



The existence of a scalar - isoscalar meson a with 7r~7r + decay was first proposed by 
Schwinger in 1957 in a field theory of strong interactions IJ. First experimental evidence for 
the existence of a meson with a mass near p° came in early 1960's from the measurements of 
forward - backward asymmetry in tt~tt + distribution in 7r~p — > 7r~7r + n production 0,[|f|,[5]|| . 
It was expected that the <r(750) resonance would show up prominently in 7t~p — > 7r°7r°n 
production where the p° does not contribute. However the measurements of this reaction 
at CERN in 1972 found no clear evidence for cr(750) @. Furthermore, in 1973, Pennington 
and Protopopescu used analyticity and unitarity constraints on partial wave amplitudes 
in 7T7T — > 7T7T scattering (Roy equations) to show that a narrow a(750) resonance cannot 
contribute to tot scattering From these facts it was concluded that a(750) does not exist 
and in 1974 Particle Data Group dropped this state from its listings. 

In 1972, van Rossum and his spin physics group at Saclay reported the first measure- 
ments of recoil nucleon polarization in ttN tcN elastic scattering made at CERN at 6 
GeV/c 0. The resulting complete set of observables enabled the first model independent 
amplitude analysis of a hadronic reaction [|10j. The results invalidated all Regge models 
with the exception of the Barger-Phillips p + p' model. The same group later measured 
at CERN the polarization in K N charge exchange and determined the tensor a<i exchange 
amplitudes from this data ||iT |. Not a single Regge model was in agreement with the mea- 



sured amplitudes. These CERN experiments established that experimental determination 
of helicity amplitudes in hadron scattering using measurements with spin is crucial for our 
understanding of hadron dynamics. 

In 1978, Lutz and Rybicki extended the concept of amplitude analysis to pion pro- 



duction processes. They showed [12| that almost complete amplitude analysis of reactions 



ttN — > 7r + 7r~iV and KN — > K + n~N is possible from measurements in a single experiment 
on a transversely polarized target. More recently it was shown [|i~3l|l4! that amplitude anal- 
yses of reactions ir~p — > 7T ir°n and ir~p — > rj-K^p are also possible from measurements on 
transversely polarized target. The work of Lutz and Rybicki opened a whole new approach 
to hadron spectroscopy and hadron dynamics by enabling us to study the production of 
resonances on the level of spin amplitudes rather than spin-averaged cross-sections. The 
measured spin amplitudes provide more than information about the resonance parameters. 
They also connect the resonances to the dynamics of hadron production. The behaviour of 
resonances on the level of spin amplitudes thus reveals new information about the production 
mechanism itself. 

The pion production on polarized targets was measured at CERN in n~p — > 7i^ii + n 
at 17.2 GeV/c [0,0,0,|n||lg] , in vr+n -> n+n-p at 5.85 and 11.85 GeV/c 0,[2T],|22| and 



in K + n — > K + ir p also at 5.98 and 11.85 GeV/c |23| , |24]j25[| . Recently measurements of 



ix~p — > 7i~7i + n on polarized target at 1.78 GeV/c were made at ITEP |26|| . Pion production 
in pp — > p-K + n using a polarized proton beam was measured at ANL at 6 GeV/c |]27| . 

The CERN measurements of n~p — > 7i^7i + n and 7i + n — > tt + tt~p on polarized targets 
reopened the question of the existence of cr(750) scalar meson. Evidence for a narrow 
cx(750) was found in amplitude analyses of n~p — > -K^-K + n at 17.2 GeV/c at low momentum 
transfers —t = 0.005 — 0.20(GeV/c) 2 and in n + n — > tx+tx~p at 5.98 and 11.85 GeV/c at larger 
momentum transfers —t = 0.2 — 0A(GeV/c) 2 p8| , |29| , |30| , |3ll . New evidence for cr(750) comes 
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from the amplitude analysis of the ITEP data at 1.78 GeV/c and -t = 0.005-0.20(Ge^/c) 2 
2|. The best fit in Q gives m a = 753 ± 19 MeV and T a = 108 ± 53 MeV. These values 



are in agreement with the ITEP results which give m a = 750 ± 4 MeV and T a = 119 ± 13 
MeV. 

For dipion masses below 1000 MeV the dipion system is produced in spin states J = 
(S— wave) and J = 1 (P— wave). There are two S— wave and six P— wave production 
amplitudes. The amplitude analysis is carried out in terms of normalized recoil nucleon 
transversity amplitudes and determines their moduli and cosines of certain relative phases 
]T2| , |22|] . The two S— wave transversity amplitudes S and S can be expressed in terms of 
nucleon helicity nonflip amplitude So and helicity flip amplitude Si. The amplitudes S and 
Si exchange ai and ir quantum numbers in the t— channel, respectively. They are of physical 
interest since the residue of the pion pole in the flip amplitude Si is related to the partial wave 
in 7T7T — ► 7i7i scattering and the residue of the ai pole in the nonflip amplitude is related to 
the partial wave in nai — > titi scattering. The helicity amplitudes cannot be determined from 
the measured transversity amplitudes since the CERN data on polarized targets do not give 
information about the relative phases between the transversity amplitudes. Experimentally, 
such information can be obtained only in measurements of recoil nucleon polarization. 

The measurements of ir~p — > 7i^7i + n on unpolarized target at CERN [[32] were used in 
model dependent determinations of titi phase shifts [p3| , |34| , p5| , |36| , p7| . Various analyses of tttc 



scattering used these CERN-Munich phase shifts to estimate the mass and the width of the 
a meson with differing results [||,||,|3,|l|,||,|3,||,||,|§] . The CERN-Munich phase shifts 



were obtained using an enabling assumption that all ai exchange nonflip amplitudes are zero 
in pion production. This assumption was invalidated by the CERN measurements of pion 
production on polarized targets which revealed large and nontrivial ai exchange amplitudes 
for any dipion spin P^|5D| . 

In 1997, Kaminski, Lesniak and Rybicki supplemented the CERN data on polarized 
target by simplifying assumptions about the absolute phases of P— wave amplitudes which 
enabled them to isolate the helicity flip amplitude Si and to estimate from it the the S— wave 
phase shift 5$ in titi scattering They found four solutions for 5q. The two "steep" 

solutions showed evidence for a narrow <r(750) state. However the obtained inelasticity 
7/0 violated unitarity and contradicted certain Air production data ^7|. The authors 



concluded that a narrow er(750) is excluded in irir scattering and thus also in pion production 
tt~p — > 7T~7T + n. Remarkably, their "flat-down" solution was similar to the old CERN-Munich 
phase shift 5®. 

In this work we present a new model independent amplitude analysis of the CERN data 
on 7t~p — > 7T~7r + n on polarized target using a Monte Carlo method. The analysis extends 
the range of dipion mass from 600 - 900 MeV used in the previous study [|29|,|3(| to 580 - 
1080 MeV. The new results for the unnormalized spin "up" amplitude \S\ 2 H show a narrow 
resonance structure below 880 MeV followed by an enhancement above 900 MeV and a 
dramatic dip ~ 1000 MeV. The spin "down" amplitude \S\ 2 H shows a broader structure 
around 700 MeV and a dip at ~1000 MeV. Here £ = d 2 a/dmdt is the integrated cross- 
section taken from Ref. |32| . 

In our next step we supplement the CERN data with an assumption of analyticity of 
pion production amplitudes in the dipion mass m 2 . The production amplitudes then satisfy 
generalized fixed s and t dispersion relations in dipion mass |48|]. In our finite mass range 
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the analyticity in dipion mass allows us to write the nucleoli transversity amplitudes as 
a sum of Breit-Wigner amplitudes for cr(750) and / (980) with complex coefficients and a 
complex coherent background. We use this parametrization to perform simultaneous fits to 
the measured amplitudes |S| 2 £ and |S| 2 £ to determine the transversity amplitudes S and S. 
We find two different fits, A and B, with the same x 2 / do} an d the same resonance parameters 
for cr(750). The solution average values are m a = 778 ± 16 MeV and r a = 142 ± 33 MeV 
for both fits. 

The fixed s and t dispersion relations for transversity production amplitudes determine 
their absolute phase. We show that the phases of the fitted amplitudes must be the absolute 
phases of the production amplitudes - up to an overall phase common to all amplitudes. The 
CERN data supplemented by the analyticity of production amplitudes in dipion mass thus 
allow a complete determination of the helicity amplitudes from the fitted transversity am- 
plitudes without the measurements of recoil nucleon polarization. This effect of analyticity 
to reduce the number of spin measurements required for amplitude determination was first 
observed in an analysis of ttN charge exchange data. When the data on polarized target for 
n~p — > 7T°n were supplemented by fixed t dispersion relations for the nonflip and flip helic- 
ity amplitudes, the fitted helicity amplitudes reproduced the measured helicity amplitudes 
obtained from measurements of recoil nucleon polarization in irN scattering | |49| |. 

There are two solutions for the helicity amplitudes corresponding to the sign ambiguity 
of the cr (770) contribution to the moduli |S| 2 E and |S| 2 £. In the "down" solution the 
contribution of cr(770) to the 7r exchange helicity flip amplitude Si is suppressed. Almost all 
contribution of <r(770) to pion production is in the a± exchange helicity nonflip amplitude So- 
la the "up" solution the situation is reversed with almost all contribution from <r(770) being 
in the 7r exchange amplitude Si. We show that unitarity in tttt scattering clearly excludes 
this "up" solution. Futhermore, we show that the "down" solution - and thus the evidence 
for the narrow cx(770) resonance - is in full agreement with unitarity in both 7r~7r + — ► 7r~7r + 
and n^al — ► 7r~7r + scattering. 

We conclude that <r(770) is suppressed in 7r~7r + — > 7r~7r + scattering in agreement with 
the CERN-Cracow phase shift analysis and the analyses of tttt scattering using Roy 
equations ||,[5(],[)1]]. The suppression of cr(770) in 7r~7r + — ► 7r~7r + scattering however does 
not imply its absence in pion production processes ir~p — ► n^7i + n and 7i + n — > 7r + 7r~p. 
Instead, <r(770) contributes significantly to pion production due to its strong presence in the 
ir~af — > 7r~7r + scattering in the nonflip amplitude So- 

There are four "down" solutions corresponding to the four combinations i,j = 1, 2 

of the two solutions for the moduli |5| 2 E(i) and |5'| 2 E(j). In the solution combinations 
(1, 1) and (2, 1) the narrow cr(770) manifests itself as a broad resonant structure at 720 
MeV in the helicity flip |Si| 2 mass spectrum and thus also in the S- wave in 7r~7r + — > 7r _ 7r + . 
The metamorphosis of a narrow <r(770) in the helicity nonflip amplitude So into a broad 
resonant structure in the helicity flip amplitude Si in these solutions is a new phenomenon 
with important connections to symmetries of QCD. 

In a related paper we show that Weinberg's mended symmetry |k| selects solutions 
(1, 1) and (2, 1). Ellis and Lanik derived a relation between the mass m a and the partial 
width r(er — > 7r~7r + ) in an effective field theory with broken scale and chiral symmetry [531 . 
In Ref. |5_2| we show that Ellis Lanik relation selects solution (1, 1) and imparts the narrow 
cx(770) resonance a dilaton-gluonium interpretation. The solution (1, 1) gives m a = 769 ±13 



4 



MeV and r a = 154 ± 22 MeV. We thus refer to the a resonance as er(770). 

The paper is organized as follows. In Section II we present and discuss the results and 
reliability of our new amplitude analysis. In Section III we introduce the analyticity of 
production amplitudes in dipion mass. In Section IV we use it for the parametrization of 
unnormalized .S'-wave transversity amplitudes in terms of <r(750) and /o(980) Breit-Wigner 
amplitudes and coherent backgrounds and present the results of simultaneous fits to ampli- 
tudes \S\ 2 Ti and (S'pS. In Section V we use the fitted transversity amplitudes to determine 
the helicity amplitudes S and Si. In Section VI we show that the" up" solution is excluded 
by the unitarity in scattering while the "down" solution is in agreement with unitarity 
in both 7T7T and nai scattering. In Section VII we comment on comparisons with tc tt° mass 
spectra. The paper closes with a summary in Section VIII. 

II. AMPLITUDE ANALYSIS 

A. Data and results for transversity amplitudes 

The high statistics CERN-Munich measurement of n~p — > 7r~7r + n at 17.2 GeV/c on 
polarized target was reported in four data sets with kinematics given in the following Table: 



Set —t m Am Ref. 

(GeV/c) 2 (MeV) (MeV) 

1" 0.005-0.20 600-900 20 [15] 

2 0.01-0.20 580-1780 40 [16] 

3 0.005-0.20 580-1600 20 [17,19] 

4 0.20-1.00 620-1500 40 [18] 



In the Table Am is the size of the mass bins. There is only one t-bin covering the whole 
interval of indicated —t. 

For invariant masses below 1000 MeV, the dipion system in reactions ttN — > 7t~tt + N is 
produced predominantly in spin states J = (S-wave) and J = 1 (P-wave). The exper- 
iments on polarized targets then yield 15 spin-density-matrix (SDM) elements, or equiva- 
lently 15 moments, describing the dipion angular distribution |l~2] , |22|] . The measured normal- 
ized observables are expressed in terms of two S-wave and six P-wave normalized nucleon 
transversity amplitudes. In our normalization 

|S| 2 + |S| 2 + \L\ 2 + \L\ 2 + \U\ 2 + |!7| 2 + |Af + |iV| 2 = 1 (2.1) 

where A = S, L, U.N and A = S, L,U, N are the normalized nucleon transversity amplitudes 
with recoil nucleon transversity "down" and "up" relative to the scattering plane. The S'- 
wave amplitudes are S and S. The P-wave amplitudes L, L have dimeson helicity A = 
while the pairs U, U and N, N are combinations of amplitudes with helicities A = ±1 and 
have opposite t-channel-exchange naturality. The unnatural exchange amplitudes L, L, U, 
U receive contributions from 'V and "ai" exchanges. The natural exchange amplitudes N, 
N are both dominated by "02" exchange. 

Amplitude analysis expresses analytically the eight normalized moduli 

|S|,|S|,|L|,|L|,|f/|,|F|,|V|,|Ar| (2.2) 
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and six cosines of relative phases 



cos(7 5L ),cos(7 5 r;),cos(7 



LU, 



(2-3) 



cos(7 5L ),cos(7 5C/ ),cos(7 



LUj 



in terms of the measured observables |12 , [15 , 22 1 . In (2.3) cos 7 and cos 7 / are cosines of relative 
phases between pairs of amplitudes with opposite transversity (e.g. between S and L, and 
between S and L). The relative phases are not independent 



lSL ~ 1SU + 1LU = 



(2.4) 



Isl ~ Isu + 1lu = 
and the cosines thus must satisfy relations 

cos(7 5L ) 2 + cos(7 S(/ ) 2 + cos(7 Ll/ ) 2 - 2 cos(7 5L ) cos(7 SC/ ) cos(7 



LU, 



= 1 



(2.5) 



cos(7 5L ) 2 + cos(7 5C/ ) 2 + cos(7 LC7 ) 2 - 2 cos(7 5L ) cos(7 5C/ ) cos{^ LU ) = 1 

There are two similar solutions in each (m, t) bin [0,^]. However in some (m, t) bins the 
solutions are unphysical: either a cosine has magnitude larger than 1 or the two solutions 
for moduli are complex conjugate with a small imaginary part. Unphysical solutions also 
complicate error analysis. Two methods are used to find physical solutions for amplitudes 
and their errors. They are x 2 minimization method and Monte Carlo method. The two 
methods are described in P^ppjrn ] . 

All amplitude analyses in [|l^,|l^,|l^,|l^,|l^] use x 2 minimization method to find solution 



for amplitudes and their errors. Only Ref. [ 15] reports the results for amplitudes for the Set 
1. In [[HJ and [I9[ the ratios |S|/|S| and S-wave intensity Is are given, so it is possible to 
reconstruct the amplitudes \S\ 2 T, and \S\ 2 T> for the Sets 2 and 3. These S- wave amplitudes 
for the Set 2 are shown in Fig. 5 of Ref. |3(J and in Fig. 5 of this paper for the Set 3. 



In Ref. |29| we used Monte Carlo method for finding the physical solutions for amplitudes 
and their errors in 7r _ p — > 7r~7r + n using the data Set 1 for dipion masses in the range 600-900 



MeV. In Ref. |30| we show that the results for the S'-wave amplitudes obtained using Monte 
Carlo method agree with the results using x 2 minimization method in [|15|| . 

In this work we report results of Monte Carlo amplitude analysis of 7i~p — > 7i^7i + n 
reaction using the data Set 3 in the mass range 580-1080 MeV in order to study the <r(750) — 
/o(980) interference. Above 1000 MeV our analysis is only approximate as we neglect the 
.D-wave. However, below 1080 MeV the -D-wave contribution is still small and its neglect 
does not affect the structure of S'-wave amplitudes. This is confirmed by direct comparison 
with S-wave amplitudes obtained in Jl9| using the x 2 method with D-wave included above 
980 MeV (see Fig. 3b below). The Monte Carlo method is based on 40,000 selections of 
spin-density matrix elements within their errors in each (m, t) bin. No physical solution was 
found in 3 bins - 650 MeV, 850 MeV and 930 MeV. 
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The amplitude analysis is carried out for normalized amplitudes \A\ 2 and \A\ 2 , A = 
S, L, U, N. However, from the spectroscopic point of view the relevant information is con- 
tained in the unnormalized amplitudes |^4| 2 £ and |A| 2 £ where £ = d 2 a /dmdt is the reaction 
cross-section. It is the unnormalized amplitudes |v4| 2 £ and |^4| 2 £ which represent the direct 
spin-dependent contributions to the mass distribution d 2 a /dmdt. The unnormalized moduli 
|A| 2 S and |y4| 2 S were calculated using the £ = d 2 o /dmdt from Fig. 12 of Ref. [f32 



The results of our new Monte Carlo amplitude analysis of ir~p — > n~n + n in the mass 
range 580-1080 MeV are shown in Fig. 1. First we observe is that the p° peak in da /dmdt 
is not uniformly reproduced in all P-wave amplitudes. In fact, the amplitudes |iV| 2 £ and 
|[/| 2 E show considerable suppression of p° production. The shapes of mass distributions 
with opposite nucleon transversities show considerable differences in both S- and P-wave 
amplitudes. The amplitudes |S| 2 £, |£| 2 S, |?7| 2 S with recoil nucleon transversity "down" 
are smaller and broader than the amplitudes |S| 2 £, |£| 2 £ and |t/| 2 S with recoil nucleon 
transversity "up". The opposite is true for the natural exchange amplitudes |iV| 2 £ and 
|iV| 2 £. 

The S-wave amplitude |S| 2 £ shows a resonant behaviour below 880 MeV in both solu- 
tions. The relative phase 7 Si between S and L amplitudes is near zero in Solution 1 and 
a small constant in Solution 2. Since |-£/| 2 £ clearly resonates, the amplitude |S| 2 £ must 
resonate as well. The amplitudes |S| 2 £ and \L\ 2 I1 are both broader and their relative phase 
7sx is again near zero in Solution 2. This strongly suggests that the amplitude |S| 2 £ also 
resonates. 

The amplitude |S| 2 £ shows an enhancement between 880 and 980 MeV. Both solutions in 
both S-wave amplitudes dip at 1010 MeV. This dip is accompanied with a dramatic change 
of phase 7 5L and ^ su above 1000 MeV. This behaviour is interpreted as evidence for a 
narrow resonance / (980). Its contribution must be taken into account in fits to amplitudes 
|S| 2 £ and |S| 2 £. 



B. 5-wave intensity 

The spin-averaged S'-wave intensity is defined as Is = i\S\ 2 + |5| 2 )£. Since there are 
two independent solutions for the amplitudes ^(i)) 2 and ^(j)! 2 , i,j = 1,2, we obtain four 
solutions for Is which we label as follows 

I s (i,j) = (\S(i)\ 2 + \S(j)\ 2 )Z (2.6) 

The results are shown in Figure 2. We see that the solutions is(l, 1) and 7^(2, 1) are clearly 
resonating at ~ 770 MeV while the solutions is(l, 2) and is(2, 2) in addition to a resonant 
structure show enhancement in the range 880-980 MeV. For comparison we show in Figure 
3 the results for Is in 7i~p — ► 7r~7r + n at 1.78 GeV/c from the measurements at ITEP also 
at —t = 0.005 - 0.20 (GeV/c) 2 . The data show a clear resonant behaviour at ~ 750 MeV 
and some energy dependence of the S'-wave intensity. The S-wave intensities in CERN 
measurements of ir + n — > 7T + 7T~p at 5.98 and 11.85 GeV/c at larger momentum transfers 
—t = 0.3-0.4 (GeV/c) 2 show similarly clear resonant behaviour at ~ 750 MeV [29,30]. 
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C. Reliability of the amplitude analysis 



In Monte Carlo amplitude analysis each physical solution obtained from a Monte Carlo 
selection of input spin density matrix elements satisfies the normalization condition (2.1) 
and the phase conditions (2.5). The physical solutions are collected and each modulus and 
cosine has its own range and distribution of values. The mean values of all moduli and 
cosines of relative phases are calculated independently from the collected physical solutions 
as the average value of the distribution. If these averages are to represent the measured 
amplitudes they also must satisfy the conditions (2.1) and (2.5). 

To test the moduli for the normalization condition (2.1) we calculated the l.h.s. of 
Eq. (2.1) for the four combinations of solutions (1, 1), (2, 1), (1, 2), (2, 2). To test the cosines 
for the phase conditions (2.5) we calculated the l.h.s. of Eqs. (2.5) for both solutions 1 and 
2 for cos(t~) and for both solutions 1 and 2 for cos(7). The results are shown in Fig. 4. 

Below 850 - 870 MeV the mean values of the moduli satisfy the normalization (2.1) 
exactly in all four combinations of solutions. Above 850 MeV there are some small deviations 
except for the combination (1, 1) which is still exact. Similarly, the mean values of the cosines 
satisfy the phase conditions exactly below 870 Mev for all solutions. Above 870 MeV there 
are small deviations for solution 2 of cos(t~) and for solution 1 for cos(7) (except one point 
at 1010 MeV). The deviations for the other solutions are largest around 1000 MeV but are 
still small. 

The observed deviations are consistent with our neglect of the .D-wave in our analysis. 
From Fig. 4 we see that the moduli are less sensitive to the neglect of the D-wave than the 
S- and P- wave interference terms (the cosines). The agreement of the moduli and cosines 
with the conditions (2.1) and (2.5) is a powerful test of the self-consistency of the data and 
the reliability of our amplitude analysis, in particular below 900 MeV where the agreement 
is exact. 

In Figures 5 and 6 we summarize the results for the 5*-wave amplitudes. Figures 5 and 
6 show the results using the y 2 minimization and the Monte Carlo methods of amplitude 
analysis for the data Set 3, respectively. The results of the two methods for the data Set 



1 are shown in Figures 1 and 2 in Ref. ]3"0 ]. The two methods are clearly consistent. The 
data Sets 1 and 3 yield similar results below 900 MeV, both providing model and solution 
independent evidence for a resonant strucure near 750 MeV. Comparing Figures 5 and 6 for 
the Set 3 we note that the values of \S\ 2 Y, below 700 MeV are lower in the analysis using 
y 2 minimization method. Above 960 MeV the analysis using y 2 method includes D-wave 
contributions and the results are similar to Monte Carlo analysis which neglects the D-wave 
below 1080 MeV. 



D. Interpretation of the measured amplitudes 

We now comment on the interpretation of the measured moduli squared which will be 
important for their parametrization. To this end we start with the data analysis. 

The reconstructed and accepted events of ir~p — > 7i^7i + n are grouped in t- and m- bins. 
In a given t- and m-bin each event has its own value of t and m. The average values of t 
and m in a given t- and m-bin are t*, m*. 
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The dipion angular distribution is described in terms of unpolarized and polarized spin 
density matrix (SDM) elements. For invariant dipion masses below 1000 MeV the experi- 
ments on polarization targets determine 15 SDM elements pi, % = 1, . . . , 15 (6 unpolarized, 
9 polarized). The SDM elements pi = pi(s,t,m) depend on energy s, momentum transfer t 
and invariant mass m. The maximum likelihood optimization in data analysis of events in 
a given (t, m) bin treats all events as if they had the same t and m. The produced SDM 
Pi are assigned to the bin average values £*, m*. Hence the result of the data analysis are 
values of SDM elements pi(t*,m*) at the particular values of t and m, namely t* and rn*. 
The fact that each event in any given (i, m) bin has its own value of t and m is reflected in 
the errors on p^ at t*, m*. 

It is a common misconception to think about the measured SDM elements (or, equiva- 
lently, moments) as average values of SDM elements over the given t- and m- bin. We can 
see that this is not the case by using the mean value theorem. We get for a bin average of 
an element pi 



11" 



At Am 



Pi(t,m)dtdm = pi(t*,m*) (2.7) 



tl 



The bin average of each p\ is thus evaluated at its own values of t*, m* which are not equal 
to £*, m*. We then conclude that the measured SDM elements are simply values of pi at the 
common points t*, m*, namely Pi(t*, m*). 

The relations between SDM elements and the moduli of amplitudes and cosines of relative 



phases given in Ref. [12,15,22] are valid at any value of s, t and m. In particular, they are 



valid at the bin averaged values t*, m*. Hence the analytical solution of these equations 



I2| , [i5| , |2"2|| determines the moduli of amplitudes and cosines of relative phases also at the bin 
averaged values t*, m*. We will use this conclusion when we discuss the parametrization of 
S'-wave production amplitudes in the Section IV. 

In the CERN analysis the data sets 1-4 described above all have one t-bin and a number 
of mass bins, k = 1,2,..., N m . The average values of t and m in the k-th (t, m) bin are 
t* k) m* k . The values of t* k are very close and can be replaced by their average value t*. 



III. ANALYTICITY OF PRODUCTION AMPLITUDES IN DIPION MASS m 2 

Our starting point is the well-known |B5| dispersion respresentation of a complex function 



f(z) with simple poles at z n , n = 1, 2, . . . , N in the complex plane z, a branch cut along a 
positive real axis from a to oo and with asymptotic property |z|/(z) — > as \z\ — > oo. We 
shall also assume that the function f(z) is a real function f(z*) = f*(z). Using Cauchy's 
integral theorem and the process of contour deformation, it can be shown |55j that 



oo 



f(z) = E -^L- + - / Imf / iX ' )dX ' (3.1) 



n=l Z ~ Z n 7T 



A remarkable feature of the proof of (3.1) is that it takes place for a fixed value of z [ p5[| . As 
the result, the dispersion relation (3.1) is also valid for moving poles for which z n = z n (z). 
In such a case the residues R n in (3.1) also depend on z, i.e. R n = R n (z). Furthermore, 
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the dispersion relation (3.1) is easily generalized to include a left-hand cut and for functions 
that are not real. In the latter case Imf(x') in (3.1) is replaced by a discontinuity function 
along the cut(s). 

The dispersion relation (3.1) can also be generalized to complex functions f{z,t) that 
depenend on a real parameter t. If for values of t in some interval < t%, ti > the conditions of 
the theorem are satisfied, then the function f(z,t) will satisfy generalized fixed-t dispersion 
relations 

* R n (z,t) 1 °f Imf(x',t)dx' 
f{ z ,t) = 2^ TT + - / ) + l.h.cut{z,t) 3.2 

^[Z-Z n (z) 7T J X'-Z 

Let s = Rez and let us assume that the complex poles z n (z) have a Breit-Wigner form 

z n = s n (s) =m 2 n - im n T(s), n = 1, . . . , N (3.3) 
We define Breit-Wigner amplitudes 

s -m z n + im n T n (s) 
and redefine the pole residues to obtain fixed-t dispersion relations 

N 

f(s, t) = J2 Rn(s, t)a n {s) + J(s, t) (3.5) 

n=l 

where R n (s,t) are the redefined pole residues and I(s,t) are the dispersion integrals over 
the the left- and right- cuts. Any necessary subtractions are included in the part I(s,t). 

An example of the parametric dispersion relations (3.5) are dispersion relations for partial 
wave amplitudes T J (s,t) in off-shell mi — > ttti scattering with off-shell mass squared t. 
Such off-shell partial wave amplitudes are part of the description of tt exchange production 
amplitudes in 7r~p — > 7r~7r + n. 

The angular distribution of the dipion state 7t~tt + in n~p — > 7r~7r + n is described by 
partial wave production amplitudes H{ Xn 0Xp (s,t,m 2 ) defined by angular expansion of the 
production amplitudes H\ nt0 \ p (s, t, m 2 , 9, (p) where the angles 6 and cf) describe the direction 



of the n~ in the ir~ir + rest frame |3Q]48|1 . Here J and A are the spin and helicity of the 
dipion state in the ^-channel; X p and A n are the s-channel helicities of the proton and the 
neutron; s is the c.m.s. energy squared, t is the four-momentum transfer squared and m 
is the dipion mass. Linear combinations of these amplitudes define nucleon helicity and 



nucleon transversity amplitudes with definite t-channel naturality [|12| , |30| . For masses below 
1000 MeV only S and P waves contribute and the participating transversity amplitudes are 
A and Z with A = S,L,U,N 110. _ 

The production amplitudes A(s,t,m 2 ) and A(s, t,m 2 ) involve kinematic singularities 
|56| in t and m 2 . At fixed s and t the relevant singularity in m 2 comes from the two-body 
phase space q = Jm 2 /A — /i 2 where /i is pion mass. Let F(s,t,m 2 ) and F(s,t,m 2 ) be the 



corresponding production amplitudes free from kinematical singularities in m 2 . Assuming 
analyticity of the amplitudes F and F in the dipion mass for experimentally accessible values 
of s and t, we can write fixed s and t dispersion relations 
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N 



F(s, t, m 2 ) = R n (s, t, m 2 )a n (s) + I(s, t, m 2 



(3.6) 



n=l 



where N is the number of contributing resonances and I(s, t, m 2 ) are the dispersion integrals. 
In a finite mass interval 4/i 2 < m 2 < m 2 M where only M resonances contribute we can write 



M 



F(s,t,m 2 ) = R n (s,t,m 2 )a n (s) + B {M \s,t, 



m 



(3.7) 



n=l 



where the complex backcground _E?( M ) 



N 



B (M) = ^ R m ( s ,t,m 2 )a m {s) + I{s,t,m 2 

m=M+l 



(3.8) 



The parametrization (3.7) of production amplitudes in terms of a sum of Breit-Wigner 
amplitudes with complex coefficients and a complex coherent background has been widely 
used in fits to production data [^7],^,^ jMj . Description of interfering resonances in nn 
scattering and in pion production n~p — > n^n + n and the related questions of unitarity 



constraints are further discussed in Ref. 48 



Fixed s and t dispersion relations in dimeson mass for production amplitudes were first 
introduced in 1973 by Hoyer and Kwiecinski |H| to relate production amplitudes in a + b — ► 
c + d + e to particle- Reggeon scattering a + R —>■ c + d. We comment on the relationship of 
analyticity in dipion mass and pion-Reggeon scattering in n~p — ► 7r~7r + n in Section VI. 



IV. DETERMINATION OF S- WAVE TRANS VERSITY AMPLITUDES FROM 
SIMULTANEOUS FITS TO |5| 2 S AND |5| 2 S 

A. Parametrization of S"-wave production amplitudes 

To understand the resonant structure of the S'-wave amplitudes |5*| 2 S and |5*| 2 S we 
performed simultaneous fits to these amplitudes using analyticity based parametrization 
(3.8) which included a <r(750) Breit-Wigner pole, / (980) Breit-Wigner pole and a coherent 
background. Following Ref. [0, we write the unnormalized S'-wave amplitudes in the form 

|S| 2 S = g|F| 2 (4.1) 



|Sf£ = q\F\ 



where q is the mass dependent phase space factor determined in [3(| to be simply the c.m.s. 
momentum of 7r~ in the 7r~7r + rest frame. The amplitudes F and F are unnormalized ampli- 
tudes which are parametrized in terms of resonance contributions and coherent background. 
For amplitude F = F(s, t, m) we write 



R a (s, t, m)a a {m) + Rf(s, t, m)af(m) + Q(s, t, m) 



(4.2) 
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where a R is the Breit-Wigner amplitude 

mRT i a i\ 

a R = — - 2 : (4.3) 

m R — m z — im R Y 

where m R is the resonant mass, R — a,f. In the following / will refer always to /o(980) 
resonance. The mass dependent width T(m) depends on spin J and has a general form 

r = r R (±)^^Ml (4.4) 

q R Dj(qr) 

In (4.4) q R = q(m = m R ) and Dj are the centrifugal barrier functions of Blatt and Weiskopf 

D (qr) = 1.0 (4.5) 

Dx(qr) = 1.0 + {qrf 
where r is the interaction radius. We recall that 

2 2 

Re a R = { mR ~ m )\a R \ 2 = w R \a R \ 2 (4.6) 
m R T 

Im a R = \a R \ 2 = H R 

In (4.2) the term Q(s, t, m) is the coherent nonresonant complex background and R a {s, t, m) 
and Rf(s, t, m) are complex coupling functions. The amplitude F has a form similar to (4.2) 
with obvious replacements R a — > R a , Rf — > Rf and Q — > Q. 

The energy variable s is fixed and will be omitted in the following. It is convenient to 
factor out the phase of R a and define 



R a = \R a \e^ = ^N s e^ (4.7) 

C = Rf/Ra = Ci + iC 2 

B = Q/R a = B X + iB 2 
Then the unnormalized amplitude (4.2) has the form 

F = ^{a a + Ca f + B}e^ (4.8) 
and the parametrization of \S\ 2 T, reads 

|S| 2 £ = qN s {[l + 2w a B 1 + 2B 2 ]H a + B\ + B\ + \C\ + C 2 ]H f + (4.9) 

+2[(w a H a + BJiwfd - C 2 ) + (H„ + B 2 )(d + w f C 2 )}H f } 
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The parametrization of \S\ 2 H has the same form as (4.9) with the obvious replacements 
N s -> N s , Bt -> 5,, d -> C i5 i = 1, 2. 

The equation (4.9) for |S| 2 £ and the analogous equation for IS'pE are valid for any s, t 
and any m< 1000 MeV. As discussed in the Section II. C, the amplitude analysis produces 
values of |S| 2 £ and |S| 2 £ at points (t*, m* k ), k = 1, 2, . . . N m . In our case t* = 0.069 (GeV/c) 2 
and iV m = 24 bins. Evaluating the r.h.s. of (4.9) and the corresponding equation for \S\ 2 H 
at (t*, ml) yields (at a given s) values 

N s (t*, ml), B t (t*, ml), d(t* , m* k ) (4.10) 



with i = 1,2 and k = 1, . . . , N m . To perform the fit to data on \S\ 2 T, and \S\ 2 H we must 
make assumptions about the mass dependence of functions Ns, B iy Ci, N s, Bi, Ci, i — 1, 2. 
We will assume that they depend only weakly on the dipion mass m and can be taken as 
mass independent in our fits. Then 

N s (t*),B t (t*),d(t*),Ns(t*),B i (t*),C i (t*) (4.11) 

are simply constants. They are the free parameters to be determined by our simultaneous 
fits to \S\ 2 Tj and \S\ 2 T, data. Notice that the parameters (4.11) are values of functions Ns, 
Bi, Ci, Ns, Bi, Ci, i — 1, 2 at t — t* and not average values of these functions over the t-bin 
as assumed in Ref. pi . 



B. Results of simultaneous fits to |5| 2 S and \S\ 2 T, 

Hadron resonance in a production process is represented as a Breit-Wigner pole in helicity 
or transversity amplitudes with the mass and width independent of helicities or transversities 
of external particles. We assume that the a resonance contributes to both transversity 
amplitudes F and F and that the difference in the shape of \S\ 2 T> and \S\ 2 H is due to the 
difference in coherent background and the cr(750) — /o(980) couplings. The measured mass 
distributions \S\ 2 H and \S\ 2 H were augmented at 650 and 850 MeV mass bins by data points 
from Monte Carlo amplitude analysis using the data Set 1 p9[ . 

We used the parametrization (4.9) for \S\ 2 H and \S\ 2 T, with common m a , Y a , mj and 
Tf. We used the programme FUMILI |k| to perform a simultaneous fit to data on \S\ 2 H 
and |£| 2 £. The two solutions for \S\ 2 T, are labeled 1 and 2, while the two solutions for \S\ 2 H 
are labeled 1 and 2. Thus there are 4 simultaneous fits (1, 1), (2, 1), (1, 2) and (2, 2). 

We used 15 different initial values of the free parameters, including fits with a fixed mass 
and width of /o(980) at mj = 980 MeV and = 40 MeV. Improved fits were obtained when 
rrif and Tf were allowed to be free parameters. Two solutions, labeled Fit A and Fit B, 
were found for the free parameters giving the same values of resonance parameters m a , Y a , 
rrif,Tf and the same values of x 2 /dof for each combination of solutions for the amplitudes. 
Most initial values produced the Fit B. 

The results of the best fits are shown in Fig. 6. The two fits A and B produce virtually 
identical curves. The fitted values of m a , F a ,mf and Tf and the associated x 2 /dof for each 
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combination of solution are given in Table I. The results for the normalization parameters 
Ns, Ns, the background parameters B iy B iy % — 1,2 and the /b(980) coupling parameters 
Ci, Ci, i = 1, 2 are given in the Tables II and III for the Fit A and Fit B, respectively. We 
observe that the free parameters for the amplitude | S'l 2 1 S are essentally the same for both 
fits A and B. For the amplitude |S| 2 £ the fits A and B differ chiefly in the values of Ns, B2 
and especially in C\ and C2. 

We see from Fig. 6 that the parametrization (4.9) reproduces well the structure of S'-wave 
amplitudes |S| 2 £ and \S\ 2 H in both solutions. The a resonance peaks are clearly visible in 
both solutions to |S| 2 £ while the presence of a resonance in |S| 2 £ is required to explain 
the broad structure observed in this amplitude. The interference of /o(980) with coherent 
background and the a resonance is responsible for the enhancements in the mass region 
880-980 MeV, the rapid decrease at 990 MeV and the dip at 1010 MeV in all solutions. 

We see in Fig. 6 that the cross-fits (2,1) and (1,2) are similar to fits (1,1) and (2,2), 
respectively. In Table I we notice that all fits have similar values of m a , T a , irtf and V f. The 
average values of mass and width of a resonance are 

m a = 778 ± 16 MeV , = 142 ± 33 MeV (4.12) 



This compares with m a = 753 ± 19 MeV and T a = 108 ± 53 MeV obtained in |I0| from 



fits only to \S\ 2 Ti without the interference with / (980). The simultaneous fit of |S| 2 £ and 
\S\ 2 Ti and the interference with / (980) thus require a slightly higher mass and a broader 
width of the a state. In the following we refer to the a meson as a(770) resonance. 

The average values of mass and width of / (980) resonance determined by our fits are 

m f = 973 ± 12 MeV , T f = 51 ± 25 MeV (4.13) 

This compares well with the Particle Data Group estimates of rrif = 980 ± 10 MeV and 
T f = 40 - 100 MeV. 

The program FUMILI calculates the x 2 contribution for each data point and x 2 /dof f° r 
each solution combination (i, j), i,j = 1, 2. The results are presented in the Table I and are 
the same for both fits A and B. The values of x 2 /dof show very good fits to all data with 
an overall average of x 2 /dof = 0.4420. 



C. The consistency of cr(770) with relative phases between S'-wave and P-wave 

amplitudes 

The moduli of S'-wave amplitudes |5*| 2 S and |S| 2 £ provide model and solution indepen- 
dent evidence for a a state below 800 MeV and we have used them to determine its resonance 
parameters. Additional evidence for the a resonance is provided by the cosines of relative 
phases between the S'-wave and P-wave amplitudes shown in Figures lc and Id. 

Below 900 MeV the cosines cos(7 5i ) and cos(7sx) are nearly equal to +1 in Solution 1 
and are close to +1 in Solution 2. The pairs of amplitudes (S, L) and (S, L) are thus in 
phase in Solution 1 and nearly in phase in Solution 2. Since the amplitudes \L\ 2 H and |£| 2 £ 
resonate near 770 MeV due to the presence of p°, a resonance with a mass and width similar 
to p° is expected in the amplitudes S and S. 
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Below 900 MeV the cosines cos(7 5(7 ) and cos(j su) are both nearly constant and close 
to -1 in both solutions. The pairs of amplitudes (S, U) and (S, U) are thus nearly 180° 
out of phase. Although the p° production is suppressed in the amplitudes \U\ 2 T, and |£/| 2 E 
compared to |£| 2 S and |£| 2 £, these amplitudes clearly resonate near 770 MeV and the 
near constant phase lag near 180° is again consistent with a p-like a state in the S-wave 
amplitudes. 

If the pairs of amplitudes (S, L) and (S, L) are in phase and the pairs (S, U) and (S, U) 
are nearly 180° out of phase, then the pairs of P-wave amplitudes (L, U) and (L, U) must 
also be nearly 180° out of phase and the their cosines should be nearly constant and close 
to -1 below 900 MeV. In Figures lc and Id we see that this is the case. From Fig. 4 we 
see that the cosines satisfy the phase conditions (2.5). As discussed in Section II. B , this 
self-consistency between the relative phases is an important test of the self-consistency of 
the experimental data and the validity of our model independent amplitude analysis. 

It is noteworthy that the relative phases in the amplitude analyses of ir~p — > 7r~7r + n at 

I. 78 GeV/c for -t = 0.005-0.20 (GeV/c) 2 @ and in ir+n -»■ n+n-p at 5.95 (GeV/c) 2 and 

II. 85 (GeV/c) 2 for —t = 0.2 -0.3 (GeV/c) 2 show the same self-consistency of phases and 
also find that the (S, L) and (S, L) are in phase in Solution 1 and nearly in phase in Solution 
2. These independent results strengthen further the evidence for the <r(770) resonance. 

This qualitative discussion relating the evidence for the a resonance with the measured 
relative phases can be translated into a quantitative analysis involving simultaneous fits to 
the moduli 

|S| 2 S, |S| 2 £, |I| 2 S, |L| 2 S (4.14) 

and interference terms 

X SL = |^||I|Scos(7 S l) = Re(SL*)Z (4.15) 



X SL = |S||L|Scos( 7 5l) = Re{SL*)Y> 

using parametrizations (4.2) for the S- wave amplitudes and similar parametrizations with 
p°(770) in the amplitudes L and L for dipion masses below 1080 MeV. The results of such 
an analysis could be used for a new determination of the S- and P- wave phase shifts 5® and 
8\ in 7T7T scattering. 



V. SEPARATION OF vr AND at EXCHANGE HELICITY AMPLITUDES. 

A. Analyticity and the phases of transversity amplitudes 5 and S 

The S- wave transversity amplitudes S and S are a linear combination of nucleon helicity 
non-flip a\ exchange amplitude So and nucleon helicity flip 7r exchange amplitude Si. The 
inverse relations are f3(| 

S = ±(S + S) (5.1) 
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In this Section the amplitudes S, S, So and Si are understood to be the unnormalized ampli- 
tudes. To determine the helicity amplitudes we need to know the phases of the transversity 
amplitudes. Experimental determination of the relative phase between S and S requires 
difficult measurements of recoil nucleon polarization [jnfl . Here we show that the analyt- 
icity of production amplitudes in dipion mass determines fully the phases of transversity 
amplitudes, up to a sign ambiguity that will be resolved using unitarity in tcti scattering. 

Fixed-t dispersion relations for helicity amplitudes in two-body scattering fix their ab- 
solute phases. Amplitude analysis of n~p — > 7r°n charge exchange using fixed-t dispersion 
relations determined the helicity helicity non-flip and flip amplitudes Fq and F\ from fits to 



differential cross-section da/dt and polarization P over a range of energies j4"5|| . We note 
that 

| = |F„|W (5.2, 
do 

P— = 2Re(F F*) 

The fitted amplitudes reproduce the measured amplitudes obtained using measurements 
of recoil nucleon polarization. This suggests that analyticity can be used to replace the 
unknown information about the recoil nucleon polarization. 

In 7i~p — > 7r~7r + n the partial wave intensity Ia and partial wave polarization Pa are 
given by equations similar to (5.2) 

I A = \A \ 2 + \A 1 \ 2 = \A\ 2 + \A\ 2 (5.3) 

P A I A = 2eRe(A Al) = \A\ 2 - \A\ 2 

where e = +1 for A = S,L,U and e = — 1 for A = N. The fixed s and t dispersion 
relations for the production amplitudes also fix their absolute phases. The analogy with ttN 
scattering suggests that using these dispersion relations to determine production amplitudes 
from polarized target data will also yield true amplitudes. 

The parametrizations (4.2) follow from the fixed s and t dispersion relations (3.7) and 
thus we expect the phases of amplitudes in (4.2) to be their absolute phases. However the 
fits are made to the moduli and the fitted transversity amplitudes Sfu and Sf# have their 
own definite phases. Let us suppose that the true amplitudes S and S differ in phase 

S = Sfne* , S = ~B fit e* (5.4) 

We note that 

Re(SS ) = Re(S fi t S*f it ) cos 5 — Im(S fuS*f it ) sin 5 (5.5) 
Re(SS*) = Re(SfitSf it ) cos 5 + Im(S fuS*f it ) sin S 
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where 5 = 0-0. Since Re(SS) = Re(SS*) we find sin 5 = which means that 

= (5.6) 

or 

= 0±7T (5.7) 

The amplitudes S, S, So and Si thus share a common phase e 1 ^. It is easy to show that this 
phase is common also to all P-wave (and higher) amplitudes A and A. The absolute phase 
e 1 ^ can thus be omitted and the helicity amplitudes S and Si can be calculated using either 
the pair (Sf it , +S fu) corresponding to (5.6) or (Sfit, —S fit) corresponding to (5.7) for the 
amplitudes S and S in (5.1). 

B. Helicity amplitudes So and S\ 

The results for the helicity amplitudes corresponding to the phases (5.6) with S = Sfu 
and S = +S fu are shown in Fig. 7. For each solution combination of amplitudes |S| 2 £ and 
|S| 2 E we present the real and imaginary parts of the helicity amplitudes So and Si and their 
moduli | S 1 2 and |Si| 2 for both fits A and B. 

The immediate and most striking observation is that er(770) is present almost entirely 
in the a\ exchange amplitude So- The contribution of cr(770) to the 7r exchange amplitude 
Si is severely suppressed. The S-wave pion production in n~p — > 7r~7r + n below 900 MeV is 
dominated by the a\ exchange amplitude S and cr(770). 

We now look at amplitudes S and Si calculated using the the pair (Sf it , —Sfit) for the 
amplitudes (S, S) in (5.1). This corresponds to the phases (5.7). In the following we briefly 
label So(+), Si(+) the helicity amplitudes obtained from the pair (Sfn,+Sfi t ) and So(— ), 
Si(— ) the helicity amplitudes obtained from the pair (Sfit, —Sfit)- From (5.1) we find 

ReS (-) = -ImS 1 (+) , ImS (-) = +ReS 1 (+) (5.8) 

PeSi(-) = +/mS (+) , /mSi(-) = -ReS (+) 

We see that the roles of S and Si have been interchanged. It is now the helicity flip 
amplitude |Si(— )| 2 which dominates the pion production and receives most of the <r(770) 
contribution. We thus encounter an analogue of the "up - down " ambiguity known from 
the CERN-Munich and CERN-Cracow determinations of hit phase shifts. In our case the 
"down" solution are the amplitudes S n (+) and the "up" solution are the amplitudes S n (— ), 
n — 1,2 in all four solution combinations (i,j),i,j = 1,2 of moduli |S| 2 S and |S| 2 S. The 
origin of the "up - down" ambiguity in our case is the sign ambiguity of S = ±S 

C. Relative contribution of cr(770) to helicity amplitudes 

To examine the quantitative contribution of <r(770) to helicity amplitudes we now write 
the transversity amplitudes in the form 
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S = q(R a a a + R f a f + Q) (5.9) 

S = q{R a a a + Rfaf + Q) 
where R a > 0, R a > for "down" solution. Here 

Ra = \Z% , Rf = R«C, Q = R a B (5.10) 

Ra = \J~N~s , Rf = RaC , Q = R a B 
are expressed in terms of fitted parameters. The helicity amplitudes then have a form 

S = q(R a0 a a + Rf a f + Q ) (5.11) 

Si = iq(R al a a + Rf\a f + Q x ) 

where the a couplings to helicity amplitudes are 

1 _ 1 — 

Rao = ~^2^^ a + ' ^ <jl = ^/2^ fJ ~~ (5-12) 

and similar equations for Rf n and Q n , n = 1,2. The results for R a o and R a \ are shown in 
Table IV. We find that the a couplings to flip amplitude S\ are small. They are largest and 
negative in the (2, 1) solution combination and smallest and postive in the (2, 2) solution 
combination. 

The smallness of R a \ suggests that <r(770) may be absent in the flip amplitude. We 
tested this assumption in modified fits A m and B m under a constraint that Ns = N s , i.e. 
with R a i = 0. The most notable difference with the fits A and B is the flat behaviour of 
| Si | 2 below 900 MeV shown in Fig. 7 for the Fit A m . The fits A m and B m are close to 
fits A and B but suffer from somewhat larger values of x 2 /d°f- This indicates that some 
contribution of <r(770) to the flip amplitude may still be necessary. 

Figure 7 shows important distinctions between the pairs of solution combinations (1, 1), 
(2, 1) and (1, 2), (2, 2) in the behaviour of |Si| 2 mass spectrum below 900 MeV. In solution 
combinations (1,1), (2,1) we observe a broad and low resonant structure centered at ~ 
720 MeV in fits A and B clearly above the nonresonating fit A m . In solution combinations 
(1, 2), (2, 2) the fits A and B show no clear resonant structure and nearly coincide with the 
nonresonating fit A m . From Table I we find that the average mass of o resonance in the 
solution combinations (1, 1), (2, 1) is ra a = 772 ± 12MeV. 

The broad resonant structure at 720 MeV seen in |Si| 2 in the solution combinations 
(1,1), (2,1) in both fits A and B indicates that a narrow resonance can contribute to a 
production amplitude without a marked resonance behaviour such as a prominent peak or a 
dip. In our case the narrow cr(770) resonance appears as a narrow resonance in the non-flip 
amplitude So but manifests itself as a broad resonant structure in the flip amplitude S\. 
This metamorphosis of a narrow <r(770) in the nonflip amplitude So into a broad resonant 
structure in the flip amplitude Si is a new phenomenon related to breaking of scale and 
chiral symmetry in QCD ||52||. 
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D. Relative contribution of /o(980) to helicity amplitudes 



Figure 7 shows also interesting distinctions between the pairs of solution combinations 
(1,1), (2,1) and (1,2), (2,2) in the behaviour of |So| 2 and |Si| 2 mass spectra above 900 
MeV. In this mass region we also find the differences between the fits A and B. 

Looking at the nonflip spectrum |So| 2 , the / (980) resonance appears as a shoulder or 
a small peak in solution combinations (1,1), (2,1) and as a pronounced narrow peak in 
solution combinations (1, 2), (2, 2). The fits A and B are distinguished by the height of the 
/o(980) peaks. 

In the flip spectrum |Si| 2 the / (980) resonance appears as a shallow dip in solution 
combinations (1, 1), (2, 1) and as a pronounced narrow dip in solution combinations (1,2), 
(2, 2) in the Fit A. In the Fit B such dips are preceeded by a peak at ~ 920 - 950 MeV 
which is more pronounced in the solution combinations (1,2) and (2,2). 

VI. CONSISTENCY OF cr(770) WITH UNITARITY IN tt--k + -> ir-ir + AND 

n~af -> ir-TT + SCATTERING 

In this Section we show that the "up" solution S n (—),n = 1,2 is excluded by the uni- 
tarity in 7T7t scattering while the "down" solution S n (+) is allowed. Consequently <r(770) is 
suppressed in 7r~7r + — > 7r~7r + reaction. 

Consider a nucleon s-channel helicity flip amplitude in 7t~p — > 7r~7r + n with a dipion spin 
J and i-channel helicity A = |30] 



A((s,t,m 2 ) = H^ _(s,t,m 2 ) (6.1) 

The pion exchange amplitudes A\ have a general form 

AJ R{(s,t,m*) + M( 
t — fi z 

where R{ is the pion pole residue and M{ is the non-pole term. The connection to tttt 
scattering is provided by a factorization hypothesis 

R{ = -FJ{s,t)TJ(m 2 ,t,0) (6.3) 

where T J (m 2 ,t,0) is the off-shell partial wave amplitude in 7r~7r + — > 7t~tt + scattering with 
the exchanged pion having mass t and spin 0. The factor Ff is positive. It contains 
the 7rNN form factor, phase space factors |3(J and angular momentum constraints. We 



have verified explicitely that the factorization condition (6.3) is consistent with unitarity 



conditions for partial wave production amplitudes in n~p — > 7i~7i + n derived in f48 |. The 
consistency requires that the helicity nonflip amplitude i?o+,o+ i n 71 V ~^ n ~P is essentially 
imaginary, positive and slowly dependent on energy variable s at small t and that it dom- 
inates the helicity flip amplitude i?o+,o-- These requirements are in agreement with data 
on Ototi^ V ^ P) through optical theorem and with measurements of target polarization 
and recoil nucleon polarization in ir~p — > ir~p from 6 to 40 GeV/c J51 |. The amplitude 
analysis of ir^p — > n^p reactions using fixed-t dispersion relations and all available spin 
measurements from 10 - 100 GeV/c is also in agreement with these requirements f65 [. 
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The off-shell amplitudes T/(m 2 ,t, 0) are related to the on-shell partial wave amplitudes 
f^{m 2 ) in 7r~7r + — »■ 7r~7r + scattering |)6| 

2?(m 2 ,t,0) = (*^) J //(m 2 ) (6.4) 

where qv(o//) and q n (on) are off-shell and on-shell pion momenta with 

4rn 2 ql(off) = (t-(m + tf){t - (m - ^) 2 ) (6.5) 

The connection of pion production to 7T7r scattering through factorization (6.3) was first 
proposed by Goebel |57J and Chew and Low [58]. The possible importance of the non-pole 



term was emphasized early by Gutay f69 |. However it is generaly accepted that the non-pole 
term is small at small t |j6| . To the extent that the non-pole term can be neglected at small 
t, the 5-wave partial wave amplitude in 7r~7r + — > 7r _ 7r + imparts a phase to the production 
amplitude S%. 

Partial wave unitarity in 7t~tt + — > 7r~7r + scattering imposes a positivity condition on the 
imaginary parts of partial wave amplitudes Imf^ > [48[] which implies positivity 



JtoT/(to 2 ,*,0) > (6.6) 

Assuming that the non-pole term is small and taking into account the positivity of the factor 
Ff, the unitarity condition (6.5) then implies the positivity condition 

ImA{{s,t,m 2 ) > (6.7) 

on the imaginary parts of the helicity flip production amplitudes Af for any J. 

In Figure 7 we see that ImS\ > for the "down" solutions for all solution 

combinations and for both fits A and B. For the "up" solution ImS\(—) = —ReSo(+). 
From Figure 7 we see that ImSi(-) is large and negative below 900 MeV for all solution 
combinations. We conclude that the" up" solution is clearly excluded by the unitarity in tttt 
scattering. 

Next we note that we can write the a\ exchange nonflip amplitudes Aq = Hq + 0+ in a 
form similar to (6.2) with a\ pole term and use a factorization similar to (6.3) to connect the 
production amplitudes Aq to off-shell partial wave amplitudes T/^m 2 , t, 1) in ir~af — ► 7r~7r + 
scattering. The positivity of ImT^ (m 2 , t, 1) from unitarity in n~a^ — > 7r~7r + then implies 
positivity ImA^ (s, t, m 2 ) > 0. This condition is in full agreement with our results for ImSo 
shown in Figure 7. This agreement indicates consistency of the "down" solution and the 
evidence for cr(770) resonance with unitarity in both 7r~7r + — > tt~tt + and 7T~af — > 7t~tt + 
scattering. 

A detailed look at ImSi in Figure 7 shows small negative values at ~ 980 MeV in the 
Fit A in two solution combinations (1,2) and (2,2). ImSo shows similar negative values 
in Fit A in all solution combinations except (2,1). This may indicate a small negative 
contribution from the non-pole terms. However we cannot exclude a possibility of a small 
violation of unitarity at this mass in the Fit A in these solution combinations. In such a 
case the unitarity would favour the Fit B. 

The production amplitudes H^ Xn 0X can be considered as helicity amplitudes in a two- 
body scattering a + b — > c(J,m) + d where the particle c has spin J and mass m. A 
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more rigorous treatment replaces the elementary particle exchange with a leading Reggeon 
exchange. The exchanged particle propagators are replaced by signature factors with a 
Regge energy dependence, and the amplitudes T/(m, t, 0) and T^ Om, t, 1) are replaced by 
pion - Reggeon scattering partial wave amplitudes T/(m, t, a n ) and (to, t, a ai ) with Regge 
trajectories a n (t) and a ai (t) fT0| , |TT]| . This formalism has the advantage that it can be 
applied also to unnatural and natural exchange amplitudes which are combinations of helicity 
amplitudes with A ^ 0. To the extent that the non-pole terms can be neglected at small 
t, the analyticity of production amplitudes can be related to the analyticity of the partial 
wave amplitudes in pion - Reggeon scattering. 

Finally we comment on the S-wave tttt phase shift 5q. In Ref. |52| we show that the 



preferred solution combination is (1, 1). We thus focus on the results in Figure 7a. Neglecting 
isospin 1 = 2 contribution to the S'-wave partial wave /° and assuming that the inelasticity 
?7q = 1 below 900 MeV, we get /° = |/q where the isospin 1 = amplitude 

/ ° = - sin 5° e< = -(sin 5° cos 5° + i sin 2 5°) (6.8) 
q q 

Neglecting further the non-pole term in (6.2), we observe from the Figure 7a that 
<5q > 90° since ReS\ < 0. In contrast, the CERN- Munich phase shift 5$ is a rising function of 
to with 5q starting well below 90° in this mass range |33|j3^j35|j36]j37|1 . Our helicity amplitudes 



are determined up to an overall phase factor e 1 ^. With <p = — | we get new amplitudes 
S n = —iS n , n = 1, 2. In this case ReSi = ImSi > and iroSi = —ReSi > 0. As the result, 
the unitarity in 7t~tt + — > 7r~7r + is still satisfied and the phase shift Sq is shifted by 90° so 
that 5$ < 90°. However now ImSo = —ReSo < below 800 MeV and the unitarity in 
n~a{ — > 7r~7r + is thus violated. All CERN-Munich determinations of tctt phase shifts work 
with amplitude Si in the s-channel ||66|| . The issue of unitarity in 7r~af —>■ 7t~tt + did not 
arise in these determinations of 7T7t phase shifts since all these analyses assume Sq = 0. 



VII. COMPARISON WITH 7r vr MASS SPECTRA 

A comparison of S-wave intensities Is from measurements of n~p — > 7T~7T + n and 7T + n — > 
7r + 7r~]9 on polarized targets at four different energies in three different experiments shows a 
clear evidence for <r(770) in all data - see Figure 2 and 3 above for 7r~p — > n^n + n and Figures 
9 and 10 of Ref. for 7r + n — > tt + tt~p. In contrast, recent measurements of n~p — > TT°7v°n 
by E852 Collaboration at BNL at 18.3 GeV/c [^| indicate a suppression of <r(770) in S-wave 
intensity in 7r°7r° production and a possible presence of cr(500) resonance. 

All these data come from very high quality measurements and the differences between 
them cannot be used to exclude one or the other. We suggest that these results are all real. 
The unexpected differences in the charged and neutral pion production data may be due to 
differences in hadron production mechanism to which resonances are closely connected and 
about which they provide relevant information. 

This view is supported by other experimental facts. The S-wave intensity in ir~p — > 
at 18.3 GeV shows significant variations in momentum transfer —t [fT2[. Similar strong 



t-dependence resembling an oscillatory behaviour was observed in the P-wave amplitudes 
in ir + n — > 7t + tt~p and K + n — > K + n~p |^TJ. The data show also important differences in 
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the .D-wave contributions to ir~p — > 7r~7r + n and 7r~p — > 7r°7r°n below 1000 MeV at small t. 
While for dipion helicity A = the intensity Ipo is near zero in both reactions, this is not the 
case for A 7^ 0. The corresponding natural and unnatural exchange D-wave intensities 
and Id- are both near zero in n~p — > 7i^n + n but make a significant nonzero contribution 
in 7T~p-> TrVn J72|. 

Additive Quark Model (AQM) provides a model of hadron scattering in nonperturbative 
region. AQM relates spin amplitudes in K + n — > K*°p and pp — > A ++ n which lead to 
relations between spin observables in K + n —>■ K + 7t~p and pp — > pii + n. These relations are 
in a remarkable agreement with experiments and demonstrate in another way the close 
relationship between resonances and the production mechanism . 

Integrated 7r°7r° mass spectra were also measured by Crystal Barrel Collaboration in 
pp — > 3tt° |73j and in pp — > 57r° and by GAMS Collaboration in central production pp — > 
Pf7r°TT°p s ||57|| . The spectra do not show narrow peaks near 750 MeV to indicate a presence 
of cr(770). However the spectra show marked differences in all these measurements 
with broad peaks at 800 MeV, 600 MeV and 500 MeV, respectively. These differences also 
indicate that production mechanism plays a significant role in the formation of the 7r 7r 
spectrum and challenge the naive expectation that all dipion mass spectra should show the 
same structure. 



Another point to consider is that scalars are probes of QCD vacuum [75]. The QCD 
vacuum participates in hadron scattering at low as well as at high energies |76 J77| and the 
differences in production of scalar mesons in production and decay processes may reveal new 
information about the structure of QCD vacuum. 



VIII. SUMMARY. 

We have performed a model independent amplitude analysis of reaction 7r~p — > 7i^n + n 
at 17.2 GeV/c extending the range of dipion mass from 600-900 MeV of the earlier analysis 
29| , |30|] to 580-1080 MeV in the present analysis. We have supplemented the CERN data by 



the assumption of analyticity in dipion mass of the production amplitudes. This allows us to 
parametrize the S'-wave transversity amplitudes S and S in terms of Breit-Wigner amplitudes 
for er(770) and /o(980) and a coherent background and to determine these amplitudes from 
simultaneous fits to the moduli |S'| 2 S and |5| 2 S. The assumed fixed s and t dispersion 
relations impart phases to production amplitudes. The phases of the fitted amplitudes 
are thus the true phases of the transversity amplitudes. This allows us to separate and 
determine the helicity nonflip and helicity flip amplitudes So and Si corresponding to a\ 
and 7r exchange, respectively. 

The sign ambiguity in the fitted amplitude S results in two solutions, "up" and "down", 
for the helicity amplitudes. The "up" solution is excluded by unitarity in -an scattering. 
The "down" solution - and thus the evidence for cr(770) resonance - is in agreement with 
unitarity in both 7r~7r + — > 7r~7r + and -K~a\ — > 7r~7r + scattering. 

In the "down" solution the <r(770) is suppressed in the 7r exchange amplitude S\ and 
thus also in 7r~7r + — > 7r~7r + scattering while it dominates the a x exchange amplitude So- 
However we find that the narrow cx(770) manifests itself as a broad resonance at ~ 720 MeV 



in the ir exchange amplitude S\ in two solution combinations (1, 1) and (2, 1). In Ref. |52 



we show how this dual manifestation of <r(770) is connected to breaking of scale and chiral 
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symmetry in QCD. The symmetries select the solution (1, 1) for which m a = 769 ± 13 MeV 
and r CT = 154 ± 22 MeV. 

The CERN measurements of pion production on polarized targets have opened a new 
approach to experimental hadron spectroscopy by making possible the study of resonance 
production on the level of spin amplitudes. The experiments also show the close connection 
between resonances and the production mechanism and may reveal new information about 
the structure of QCD vacuum. Our results emphasize the need for a dedicated and sys- 
tematic study of various production processes on the level of spin amplitudes measured in 
experiments with polarized targets. Such "amplitude spectroscopy" will be feasible at high 
intensity hadron facilities |78| , |79|| . The first high intensity hadron facility will be the Japan 
Hadron Facility at KEK. It will become operational in 2007 and the formation of its physics 
program is now in progress [BQ]. 
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FIGURES 



FIG. 1. Mass dependence of physical solutions for unnormalized moduli squared of S- wave 
and P-wave nucleon transversity amplitudes and cosines of their relative phases in reaction 
ir~p — > 7r - 7r + n at 17.2 GeV/c and momentum transfers —t = 0.005 — 0.2 (GeV/c) 2 . The re- 
sults are in the t-channel dipion helicity frame. 

FIG. 2. Four solutions for the S-wave partial-wave intensity Is in the reaction ir~p — ► 7r~7r + n 
at 17.2 GeV/c and -t = 0.005 - 0.20 (GeV/c) 2 . 

FIG. 3. Four solutions for the S-wave partial-wave intensity Is in the reaction ir~p — ► 7r~7r + n 
at 1.78 GeV/c and -t = 0.005 - 0.20 (GeV/c) 2 . Data from Ref. [26]. 

FIG. 4. Consistency of measured amplitudes with the normalization condition (2.1) and phase 
conditions (2.5). 

FIG. 5. S-wave transversity amplitudes |S| 2 X and |S| 2 £ measured in ir~p^ — > 7r~7r + n at 17.2 
GeV/c and —t = 0.005 — 0.20 (GeV/c) 2 using the x 2 minimization method in amplitude analysis 
of Data Set 3. Data from Ref. [19]. 

FIG. 6. S'-wave transversity amplitudes jSpS and |S*| 2 S measured in ir~p^ — ► 7r _ 7r + n at 17.2 
GeV/c and —t = 0.005—0.20 (GeV/c) 2 using the Monte Carlo method in amplitude analysis of Data 
Set 3 (this paper). The curves show the best simultaneous fits to the measured amplitudes |S*| 2 S 
and |5| 2 S in the four solution combinations (1, 1), (2, 1), (1, 2) and (2, 2) using the parametrization 
(4.9). The fits A and B yield virtually identical curves. 

FIG. 7. Helicity nonflip amplitude So an d helicity flip amplitude Si for Fits A and B in the 
four solution combinations (1,1), (2,1), (1,2) and (2,2). The Fit A m corresponds to the total 
absence of a (770) in the helicity flip amplitude Si. 
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TABLES 



TABLE I. Mass and width of a(770) and / (980) from simultaneous fits to |S| 2 E and \S\ 2 T,. 
The combination of solutions of |S*| 2 S and jS'pS is given in the parentheses on the left with the 
notation 1 and 2 for solutions of amplitude |5| 2 S and 1 and 2 for solutions of amplitude IS 1 ] 2 !!. 
Fits A and B give the same results. 



Fit 




L(j 


m f 


r/ 


X 2 /dof 




(MeV) 


(MeV) 


(MeV) 


(MeV) 





(1,T) 769 ± 13 154 ± 22 979 ± 12 60 ± 23 0.5395 

(2.1) 774 ± 12 121 ± 23 956 ± 16 74 ± 32 0.3423 

(1.2) 787 ±19 165 ± 45 982 ± 6 27 ± 26 0.5027 
(2,2) 780 ± 18 126 ± 40 975 ± 14 44 ± 20 0.3836 



TABLE II. Normalization, background and /o(980) coupling parameters for amplitudes \S\ 2 H 
and |5| 2 S from simultaneous fits: Results for Fit A. The combinations of solutions of |S| 2 S and 
|S| 2 £ as in Table I. 



Fit 


(IJ) 


(2,1) 


(1,2) 


(2,2) 


N s 


3.652 


3.239 


1.564 


1.478 


N s 


1.966 


1.560 


2.712 


1.826 


Bi 


0.522 


0.693 


0.997 


1.128 


B 2 


0.260 


0.183 


0.881 


0.893 


Bi 


-0.319 


-0.436 


-0.118 


-0.370 


B 2 


0.332 


0.656 


0.211 


0.633 


Ci 


-0.399 


0.033 


-0.627 


-0.782 


c 2 


0.499 


0.929 


1.961 


2.392 


Ci 


-0.796 


-0.439 


-0.838 


-0.533 


c 2 


0.065 


0.622 


0.161 


0.466 
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TABLE III. Normalization, background and /o(980) coupling parameters for amplitudes |S*| 2 S 
and \S\ 2 T, from simultaneous fits: Results for Fit B. The combinations of solutions of \S\ 2 H and 
|S| 2 £ as in Table I. 



Fit 


(1,1) 


(2,T) 


(1,2) 


(2,2) 


O 


3.484 


2.985 


1.524 


1.382 


Nq 


1.966 


1.556 


2.705 


1.829 


Bi 


0.518 


0.693 


0.987 


1.068 


B 2 


0.297 


0.280 


0.919 


1.035 


Bi 


-0.319 


-0.437 


-0.119 


-0.367 


B 2 


0.332 


0.657 


0.213 


0.632 


Ci 


-0.190 


0.244 


-0.113 


0.359 


c 2 


0.386 


0.595 


1.644 


1.424 


Ci 


-0.796 


-0.440 


-0.832 


-0.532 


c 2 


0.064 


0.624 


0.162 


0.468 



TABLE IV. Couplings of <r(770) in helicity non-flip and flip amplitudes So and Si in the "down" 
solution for Fits A and B. Solution combinations as in Table I. 



Fit 




RaO 


Ral 


Ral 


X 2 /dof 




Fit A 


Fit B 


Fit A 


Fit B 





(1,T) 2.343 2.312 -0.360 -0.329 0.5395 

(2.1) 2.156 2.104 -0.390 -0.340 0.3423 

(1.2) 2.049 2.037 0.280 0.290 0.5027 
(2,2) 1.815 1.788 0.096 0.125 0.3836 
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